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ABSTRACT 


Measurements were taken of the amount 
of soil eroded from small plots under the im- 
pact of a fixed amount of simulated rain. Un- 
der these conditions, erosion is more closely 
related to amount of cover than to any other 
site characteristic. However, the relation be- 
tween erosion and cover is strongly influenced 
by slope gradient. Regression analyses indi- 
cated that erosion is about the same on a 5- 
percent slope with 40-percent cover as it is on 
a 35-percent slope with 80-percent cover. Or- 
ganic matter is the most important soil param- 
eter affecting erodibility, but the direction 
and magnitude of its effects depend on soil 
texture. Organic matter decreases erosion of 
clay soils, but tends to increase erosion of 
sandy soils. 


ii 


INTRODUCTION 


During the past decade, the author has stud- 
ied infiltration and erosion potentials of seven 
diverse summer ranges and has related these to 
cover, soil, and topographic parameters by 
means of multiple regression analyses. Accepta- 
ble regression equations relating infiltration to 
site factors were developed for three of the 
areas (Meeuwig 1969, 1970a), but the regres- 
sion equations developed for the other four 
lacked precision and contained anomalous 
relations. 

Satisfactory regression equations relating 
erosion to site factors were developed for each 
of the seven areas (Meeuwig 1970b). Unlike the 
infiltration equations, the general relations de- 


All the study areas have herbaceous cover 
and all are grazed by livestock during the sum- 
mer, except the Davis County Experimental 
Watershed from which grazing has been pro- 
hibited for more than 30 years. Their locations 
are shown in figure 1. Areas studied were: 

1.— Great Basin Experimental Range, 
Manti-LaSal National Forest, central Utah. Ele- 
vations of the 162 study plots varied from 
7,000 to 10,000 feet. This sheep range has a 
wide variety of grass and forb species. Soils are 
mostly silty clay loams and clay loams derived 


fined by these equations did not differ greatly 
from area to area. In essence, these equations 
indicate that erosion potential (as measured in 
these studies) depends chiefly on cover and 
slope gradient; and that differences in the 
cover-slope-erosion relations, both within and 
among study areas, are attributable to varia- 
tions in soil properties, notably organic matter 
content and texture of surface soils. 

This paper reports the results of combining 
the data from all seven study areas to determine 
the general influences of cover, slope, soil tex- 
ture, and organic matter on soil stability over 
the range of conditions encompassed by these 
studies. 


STUDY AREAS 


from limestone, shale, and sandstone. 

2. — Davis County Experimental Water- 
shed, Wasatch National Forest, northern Utah. 
Elevations of the 80 study plots varied between 
8,000 and 9,000 feet. This area was the source 
of serious floods during the period 1923 
through 1930. Much of it was contour trenched 
and seeded with grass from 1933 through 1936. 
Grazing has been prohibited since 1933. Soils 
are mostly silt loams and loams derived from 
gneiss, schist, conglomerate, sandstone, and 
shale. 


nn 
& 
= 
= 
°o 
AS 
iy] 
oO 
7) 


Figure 1. — Locations of study areas. 


3. — Vigilante Experimental Range and 
Monument Ridge in the Gravelly Range, Bea- 
verhead National Forest, southwestern Mon- 
tana. Elevations of the 84 study plots varied 
between 7,000 and 9,500 feet. This is cattle 
and sheep range dominated in many parts by 
Idaho fescue (Festuca idahoensis) and in others 
by native forbs or seeded grasses (mainly 
Agropyron desertorum). Soils are mostly silt 
loams and silty clay loams derived from red 
shale, siltstone-shale, and glacial till. 

4.— Diamond Mountain Cattle Allotment 
near Flaming Gorge, Ashley National Forest, 
eastern Utah. This is an experimental grazing 
area where portions of the native sagebrush- 
grass vegetation has been replaced by intro- 
duced grass species. Soils are loams and sandy 
loams derived from sandstone. Elevations of 
the 34 plots were about 8,000 feet. 


Uniform procedures were followed on all 
460 plots involved in this analysis. The plots 
were 20 inches wide and 30.5 inches long; ap- 
proximately 1/10 milacre in area. Dortignac’s 
(1951) rain simulator was used to apply 2.5 
inches of water to these plots at a constant in- 
tensity of 5 inches per hour for 30 minutes. The 
raindrops produced by this simulator tend to 
be larger than those of actual thundershowers, 
but have lower impact velocities. All runoff 
from the plots was collected and measured. All 
eroded mineral and organic material, including 
that deposited in the collecting trough at the 
bottom of the plot frame and that suspended in 
the runoff, was ovendried and weighed. 

Density of cover on each plot was measured 
with a point frame (Levy and Madden 1933). 
First strikes at 100 evenly-spaced points were 


5. — The basalt study area north of Seven 
Devils, Nezperce National Forest, central 
Idaho. Forty-five study plots were located in 
grassy openings in open ponderosa pine stands 
at about 5,000 feet elevation. Soils are loams 
and silt loams derived from basalt. 

6.— Coolwater Ridge, Nezperce National 
Forest, central Idaho. Elevation of the 15 plots 
was about 6,000 feet. Vegetation on this deteri- 
orated subalpine range is predominantly forbs 
of low value insofar as palatability and protec- 
tion are concerned. The granitic soils are sandy 
loams and loams. 

7. — Trinity Mountains, Boise National For- 
est, southern Idaho. Forty study plots were lo- 
cated in large and small openings in coniferous 
forest at elevations of about 7,000 feet. The 
granitic soils, typical of much of the Idaho 
Batholith, are sandy loams and loamy sands. 


MEASUREMENTS 


recorded as plants by species, litter, stone, or 
bare soil. A day or two after the simulated rain 
test, vegetation was clipped at the soil surface 
and litter removed. Vegetation and litter were 
air-dried and weighed. 

The following soil parameters were among 
those measured: 


1.— Antecedent moisture content of the 
surface 2 inches of soil; 


2.— bulk density and capillary porosity of 
the surface 6 inches of soil; 


3.— organic matter content of the surface 
inch of soil by the dichromate method 
(Peech 1947); 


4.— particle-size distribution (Bouyoucos 
1962) and aggregation (Middleton 
1930) in the surface inch of soil. 
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Figure 2. — Calculated and measured erosion of the 460 plots. Study areas are identified to show 
variation associated with each area. 


Variables were subjected to preliminary 
screening during the development of the seven 
regression equations for the individual study 
areas. Several site variables that were expected 
to have a definite influence on erodibility (e.g., 
aggregation and antecedent soil moisture con- 
tent) did not appear in any of these equations, 
and so were not included in the present anal- 
ysis. Although stone cover materially affected 
erosion on some of the study areas, it did not 
prove to be important when the combined data 
were screened. 

Final screening of the combined data was 
limited to the following parameters, their inter- 
actions, and transforms: 

A.— Proportion of the soil surface covered 
by plants and litter, as determined by 
first strikes of the point analyzer. It is 
equivalent to the proportion of soil 
surface protected from direct impact 
by raindrops. 

L. — Air-dry weight of litter in pounds per 
milacre. 

G. — Slope gradient of plot in percent. 

C. — Proportion of the surface inch of soil 
composed of clay. 

D.— Proportion of the surface inch of soil 
composed of sand. 

M. — Proportion of the surface inch of soil 
composed of organic matter. 


Standard Correlation 
Variable Mean deviation with Y 
A 0.687 0.247 —=O02155 
L 1.87 222 — .533 
G 18.4 8.5 291 
Cc .240 .108 .063 
D .295 .210 .059 
M 075 .038 mee of 3) 
¥. —209 .803 -- 


The dependent variable Y is the common 
logarithm of the weight (pounds per milacre) of 
material eroded from the plots during the 


ANALYSIS 


30-minute simulated rainstorm. Its average of 
—0.239 is equivalent to a geometric mean of 
0.577 pound per milacre. The logarithmic 
transformation was used because the distribu- 
tion of erosion values was skewed. The trans- 
formation resulted in a more nearly normal dis- 
tribution and more nearly homogeneous 
variance. 

Taking into consideration the curvilinear 
and interactive nature of the relations among 
the parameters as observed during preliminary 
analyses, the following regression model was 
assumed: 


¥ = Bo + DBiXj 

in which the Xj were: A, A”, A®, A*, A®, A®, 

L, L?,, AL, A?L,. A°L, G, G’, AG, A?G, A®G, 
C, C?, D, D*?, M, M?, CD, CM, DM, M/C, MD/C, 
MC?, M2C, MD?, M2D. 


These components were screened in a com- 
puterized regression analysis designed to select 
those that contributed materially to the regres- 
sion model. The following equation resulted: 


¥ = — .6935—6.456 A?+17.483 AS 
—12.403 A°—.0582 A?L+.0306 G 
— .0217 A°?G+8.21C—10.59C? 
— 8.45M+.651 M/C—1.38CD+ 35.48 M?D. 


The equation explains 74 percent of the var- 
iance of the log of erosion. The standard error 
of estimate (0.42) is difficult to interpret be- 
cause it is logarithmic. A clearer picture of the 
deviations from regression is derived when ac- 
tual erosion and calculated erosion are plotted 
on logarithmic scales (fig. 2). 

This empirical equation should not be ap- 
plied indiscriminately to specific situations be- 
cause, after all, it is derived from measurements 
of erosion caused by a fixed amount of simu- 
lated rain on small plots in a few selected areas. 
In spite of its limitations, this equation pro- 
vides some indications as to the combined ef- 
fects of cover, slope, and basic soil properties 
on soil stability. 


COVER AND SLOPE 


nds Per Milacre) 


CALCULATED EROSION (Pou 


As found in other studies of a similar nature, 
erosion on the seven study areas is influenced 
more by cover than by any other site variable. 
The magnitude of the effect of cover on erosion 
depends on slope gradient. The curves in figure 
3 are derived from the above regression equa- 
tion; sand, clay, and organic matter contents of 
the surface inch of soil are fixed at their respec- 
tive averages of 30, 24, and 8 percent. Unfor- 
tunately, it was necessary to ignore the litter 
weight term in calculating these curves, but the 


35% Slope 


20 30 40 50 60 70 80 90 100 
PLANT AND LITTER COVER (%) 


Figure 3. — Erosion calculated as a function 
of slope gradient and percentage of the 
soil surface protected by plants and lit- 
ter. Sand, clay, and organic matter are 
fixed at their respective averages of 30, 
24, and 8 percent. 


effects of litter weight will be quantified 
shortly. 

The effects of cover are greater on steeper 
slopes. Similarly, slope becomes increasingly 
important as cover decreases. At less than 50- 
percent cover, erosion rates double for each 
10-percent increase in slope. The amount of 
cover needed to hold erosion within some spec- 
ified limit varies widely with variation in slope. 
For example, 40-percent cover on a 5-percent 
slope apparently is as effective as 80-percent 
cover on a 35-percent slope because calculated 
erosion is about 1 pound per milacre in both 
cases. 

The regression term for litter weight 
(—0.058A?L) indicates that weight of litter is 
not important unless cover is virtually com- 
plete (table 1). At 60-percent cover, the pres- 
ence of 1 pound of litter per milacre reduces 
estimated erosion to 97 percent of the value 
shown in figure 3. At 100-percent cover, the 
same amount of litter reduces erosion to 87 
percent. Twelve pounds of litter per milacre 
with 100-percent coverage reduces estimated 
erosion to one-fifth that shown in figure 3. 
Thus, litter weight serves mainly to explain var- 
iations in erosion when cover is complete or 
nearly so. So far as soil stabilization is con- 
cerned, however, ground coverage provided by 
litter is much more important than its weight. 


Table 1. — Correction factors for litter weight 


Air-dry weight of litter 


Cover ; 
(percent) Pounds per milacre 
1 2 4 8 12 
- - - Correction factors - - - 
GO) 0:97 950294 
70 296 9 0:83 
80 .94 87 .76 0.58 
90 sg) 82 .68 46 0.31 
100 .87 -76 .58 34 .20 


SOIL TEXTURE AND 
ORGANIC MATTER CONTENT 


Sand, clay, and organic matter contents of 
the surface inch of soil were arbitrarily fixed at 
their average values to calculate the curves in 
figure 3. Variation within these parameters can 
cause estimated erosion to vary from less than 
half to more than twice the values shown in 
figure 3. Effects on erodibility are produced by 
these soil components operating together in a 
complex manner only approximated by the re- 
gression model presented in this paper. How- 
ever, some general trends can be detected (see 
fig. 4) and are listed below: 

1.— If organic matter is low, clay is more 

erodible than sand. 

2. — Erodibility of clay is inversely related 
to organic matter content. 

3. — Erodibility of sand tends to increase 
with increasing organic matter, es- 
pecially if the soil contains less than 10- 
percent clay. 

4. — If organic matter content is high, sand 
can be more erodible than clay. 

Of particular interest is figure 4A which 
shows relative erodibility as a function of or- 
ganic matter and sand contents at a constant 
clay content of 10 percent. Erodibility is about 
half that of average soil if organic matter does 
not exceed 8 percent or if sand does not exceed 
30 percent. When both sand and organic matter 
exceed these percentages, a further increase in 
either sand or organic matter results in 
increased erodibility. 

The apparent adverse effects that organic 
matter has on stability of sandy soils are be- 
lieved to be real because comparable relations 
appeared in the three individual regression 
equations for Diamond Mountain, basalt, and 
Trinity Mountains study areas. This phenom- 
enon is believed to be linked with water repel- 


lence, which has been observed to occur mostly 
in sandy soils and to be related to organic mat- 
ter (Krammes and DeBano 1965). Organic 
coatings profoundly affect the interfacial en- 
ergy of soil particles. Organic coatings on soil 
particles reduce the affinity those particles 
have for water to such a point that positive 
pressure may be required to force water into 
the soil (Fink 1970). These coatings may also 
reduce the forces of attraction between sand 
particles to such an extent that they repel one 
another and increase their erodibility. 

Clay tends to reduce or eliminate the ad- 
verse effects of organic matter on sand. One or 
more of the following mechanisms may be 
responsible: 


1. — Organic matter is known to favor aggre- 
gation of clay. Perhaps the stabilizing 
effects of organic matter on clay com- 
pensate for its deleterious effects on 
sand. 


2.— Clay reduces the amount of organic 
matter available for binding to sand 
particles because it has a much greater 
surface area per unit weight than sand 
and a greater capacity to bind organic 
matter. 


3.— Clay may form the necessary link for 
the aggregation of sand particles. Al- 
though the organic molecules on ad- 
jacent sand particles form no bonds 
with each other, they may form weak 
bonds with adjacent clay particles. 


The foregoing is, of course, pure speculation 
but it points out the need for further research 
into these phenomena and suggests the possi- 
bility of the use of soil amendments to stabilize 
organic sands. 
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Figure 4. — Relative erodibility as influenced by variations in sand, clay, and organic matter content 
of the surface inch of soil. Relative erodibility is the ratio of calculated erosion corrected for 
sand, clay, and organic matter to calculated erosion at average sand, clay, and organic matter. 


CONCLUSIONS 


Unless plant and litter cover is virtually 
complete, erosion potential is strongly affected 
by slope gradient and the amount of cover re- 
quired to limit potential erosion to some speci- 
fied amount is substantially greater on steep 
slopes than on gentle slopes. 

Soils vary in their susceptibility to erosion. 
The most erodible soils are those with high 
clay, low sand, and low organic matter con- 
tents. The combinations that appear to be least 
erodible are soils that have: 

1. — High clay, low sand, and high organic 

matter contents; 

2. — low clay, high sand, and low organic 

matter contents; or 

3. — low clay and low sand contents, regard- 

less of organic matter content. 

Apparently, the generalization that organic 
matter favors soil stability is not universally 
true. Erodibility may increase with increasing 
organic matter content in very sandy soils. This 
poses a difficult problem for management of 
such soils; sufficient cover must be maintained 
but at the same time, excessive accumulation of 
soil organic matter discouraged. 
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a Boise, Idaho 
7 oe Bozeman, Montana (in cooperation with 
Montana State University) 
Logan, Utah (in cooperation with Utah State 
University) 
Missoula, Montana (in cooperation with Uni- 
; versity of Montana) 
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=e University of Idaho) 
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